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ABSTRACT 

A large drift in the rotation rate of Titan observed by Cassini provided the 
first evidence of a subsurface ocean isolating the massive core from the icy crust. 
Seasonal exchange of angular momentum between the surface and atmosphere ac- 
counts for the magnitude of the effect, but observations lag the expected signal by 
a few years. We argue this time lag is due to the presence of an active methane 
weather cycle in the atmosphere. An analytic model of the seasonal cycle of 
atmospheric angular momentum is developed and compared to time-dependent 
simulations of Titan's atmosphere with and without methane thermodynamics. 
The disappearance of clouds at the summer pole suggests the drift rate has al- 
ready switched direction, signaling the change in season from solstice to equinox. 

Subject headings: planets and satellites: individual (Titan) 



1. Introduction 

Models of the interior of Titan for plau- 
sible ammonia concentrations predict the 
existence of a water ocean under a thin 
i ce I crust of pe rhaps 100 km thickness 
(jSohl et al.l 120031 ) . A subsurface ocean 
would dynamically isolate the surface from 
the more massive core, effectively lower- 
ing the moment of inertia of the surface. 
For the expected magnitude of seasonal 
angular momentum change in the atmo- 
sphere, a thin ice crust experiences mea- 



surable changes in length-of-day, or equiv- 
alently spin rate, in a matte r of a few years 
flTokano fc NeubauerllJOOJ ). 

Cassini RADAR spin rate measure- 



ments (jStiles et al.l 120081 ) , seasonal winds 
fr om an atmospheric mo del prediction 



(ITokano fc Neubauerl 120051) and Hu ygens 



wind observations (IBird et al.l 120051 ) were 



used to infer the mo ment-of-inertia (MOI) 



of Titan's surface (ILorenz et al.l 120081 ). 
The surface MOI is well below that al- 
lowed for any solid internal configuration, 
but is consistent with models of Titan's in- 



terior that with a thin su rface shell oyerlay - 
ing a subsurface ocean (jSohl et al.l 120031 1 . 
Seasonally changing zonal winds in the 
atmosphere torque the freewheeling sur- 
face, introducing a forced libration about 
synchronous rotation. The forced libra- 
tion of the surface, however, unexpectedly 
lags seasonal forcing by a few years. Here 
we suggest this time lag signals an active 
methane weather cycle in the atmosphere. 
Titan's atmosphere is known to super- 
rotate, with angular momentum exceed- 
ing that of solid body rotation. Strato- 
spheric zonal winds exce ed 100 m/s in 
the prograde direction ([Hubbard et al. 



199.4 lAchterberg et al.l l2008h . Between 



10 and 60 km, the troposphere weakly 
super rotates with win ds limited to 10-20 
m/s flBirdet al.l l2005h . The lowest ~10 
km at low latitudes show slight retro- 
grade motion, which is consistent with a 
seasonal, cross-equatorial circula tion, or 
Hadl ey cell, driven by convection (jTokano 
20071 ). Winds in the lowest kilometer tur n 
prograde again (IKarkoschka et al.l 120071 ). 
which is corroborated by the directionality 



of eq uatorial dune features (ILorenz et al. 



20061 ) J*l Atmospheric models suggest angu- 
lar momentum is seasonally redistributed 
in the atmosphere by a cross-equatorial cir- 
culati on that change s direc t ion with sea- 
sons ( 



1999; 



Hourdin et al. 19951; Tokano et al. 



Tokanol 120051: iRannou et all 120061: 



Mitchell et al.ll2006l : lTokandl2007l : Mitchell 
2008h. 

Though stratospheric zonal winds dwarf 
those in the troposphere, very little mass is 



*The presence of prograde winds at the equatorial 
surface is mysterious given frictional coupling with 
the surface. Explanation of this phenomenon is 
beyond the scope of this work. 



being advected in the stratosphere. Zonal 
winds in more than 90% of the atmospheric 
mass are below 20 m/s. The lowest one- 
third of the atmospheric mass is subject 
to the largest seasonal change in angular 
momentum due to frictional torque from 
the surface acting on a seasonally reversing 
Hadley cell. In the following, we develop 
a theory for the magnitude of seasonal 
change in atmospheric angular momentum 
(AAM) and its effect on the spin rate. We 
then develop a full time-dependent model 
of this process and show by comparison 
with numerical simulations that a seasonal 
time lag of a few years is a natural conse- 
quence of methane thermodynamics in Ti- 
tan's atmosphere. We conclude by offering 
predictions for current and future Cassini 
observations. 

2. A model for seasonal atmosphere- 
surface angular momentum ex- 
change 

Here and throughout we assume ax- 
isymmetry of the atmosphere; by doing 
so we preclude the development of super- 
rotation. Asymmetric surface heating by 
sunlight drives a convective updraft at a 
certain latitude. Surface friction imparts 
angular momentum to converging, surface- 
level air, which is then lofted into the upper 
troposphere by the convective updraft. At 
the top of the convective layer, air diverges 
along meridians while conserving angular 
momentum. The Hadley cell homogenizes 
the angular momentum in this upper, di- 
vergent branch of the flow. The triple re- 
quirements of angular momentum and en- 
ergy conservation and cyclostrophic bal- 
ance (jVallisI I2006J ) allow prediction of the 
width and strength of the Hadley circula- 



tion (IHeld fc Houl Il980l ); this model pre- 
dicts a global Hadley cell on Titan. Air 
sinks at the poles and returns along the 
surface where friction resets the angular 
momentum to the surface value. 

2.1. Atmospheric angular momen- 
tum 

In our model, the maximum achievable 
specific angular momentum at any location 
in the Hadley domain, i.e. the convective 
layer of the atmosphere, is that of the equa- 
torial surface 



l max = ^ 



T 



where Q is the rotation rate and a is the 
planetary radius. The angular momen- 
tum of the up draft is nearly zero when the 
convective updraft is located at the poles. 
Provided the updraft reaches the poles at 
or near solstices, the maximum change in 
specific angular momentum is equal to the 
maximum, Equation [TJ If we further as- 
sume the Hadley cell distributes this angu- 
lar momentum evenly over a layer of thick- 
ness Ap a , the maximum change in atmo- 
spheric angular momentum is, 



AL 



at in 



M n na 2 



(2) 



where 

M a = 47ra 2 A Pa /g (3) 

is the mass of the layer where the Hadley 
cell homogenizes angular momentum (as- 
suming hydrostasy). 

2.2. Effect on spin rate 

Conservation of angular momentum 
dictates that the surface spin rate must 
change in response to changes in atmo- 
spheric angular momentum. Seasonal 



changes in atmospheric angular momen- 
tum can measurably alter the spin rate of 
the surface, particularly if there exists a 
subsurface ocean iso lating the icy litho- 
sphere from the core (jTokano fc Neubauer 



2005 



Lorenz et al.l 120081 ). The moment of 



inertia for a thin, isolated shell of thickness 
Ar; is 

h ^ \M % a 2 (4) 



where 



Mi ~ Aira 2 piAri 



(5) 



is the mass of the shell and pi is the density 
of ice I. The maximum seasonal change in 
the spin rate of this shell is 



A tt r 



~AL atm / I s 
2 Mi 



(6) 



2.3. Maximum angular momentum 
exchange and spin rate drift 

The convective layer on Titan extends 
from the surface (1500 mbar) to roughly 
1000 mbar. If the upper branch of the 
Hadley cell extends through the top half of 
this layer, Ap a = 250 mbar; taking g = 1.35 
m/s 2 , a = 2575 km, and fi = 4.5 x 10~ 6 , 
the maximum seasonal change in atmo- 
spheric angular momentum in the Hadley 
domain is AL atm ~ 2.5 x 10 25 kg m 2 /s. 
This value is comparable to the change pre- 
dicted in a full general circulation model 
(GCM) dTokano fe Neubauer! IJOOsh . sug- 
gesting the seasonal cycle of the Hadley 
cell is dominating changes in the atmo- 
spheric angular momentum. If we assume 
Pi = 917 kg/m 2 , Ari = 70 km and Ap a 
= 150 mbar, then M a /M { ~ 4 x 10 -4 
and the maximum seasonal spin rate drift 
is AQ max ~ 2 deg/year, which is again in 



rough c orrespondence with the GCM pre- 
diction (ILorenz et al.ll2008l ). 



2.4. Time-dependent angular mo- 
mentum exchange 

The arguments presented so far can be 
generalized into a time-dependent predic- 
tion of spin rate change for comparison 
with observations. We make the follow- 
ing assumptions: 1) the convective updraft 
of the Hadley cell follows in-phase with 
the latitude of maximum diurnal-average 
solar forcing, ip 1 and 2) the Hadley cell 
instantaneously resets atmospheric angu- 
lar momentum to the value of the sur- 
face at if\. With these assumptions, the 
time-dependent atmospheric angular mo- 
mentum is 



L atm ((pi(t)) = Const.+-M a Qa 2 cos 2 ((pi(t)) 

(7) 
where the constant results from the plan- 
etary and stratospheric angular momenta 
that we assume change very little. All of 
the time dependence enters through the 
latitude of the updraft, ipi. 

Figure [1] shows the time-dependent 
model prediction of atmospheric angular 
momentum for previously chosen parame- 
ters with a constant subtracted (solid line). 
The latitude of maximum diurnally aver- 
aged solar forcing discontinuously switches 
from mid summer latitudes to the summer 
pole near solstices, during which time the 
model atmospheric angular momentum is 
a minimum. Figure [2] shows the model pre- 
diction of spin rate drift (with the average 
subtracted) assuming the surface is a 70 
km thick ice I shell (solid line). The bold 
l ine shows the m easured spin rate drift 
(jStiles et al.ll2008l ). The model clearly pre- 



dicts the correct magnitude, but fails to 
produce the observed time lag. 

3. The role of methane thermody- 
namics in producing a time lag 

The updraft of the Hadley cell is lo- 
cated where convection is favorable, which 
on Earth is nearly steadily located at warm 
equatorial oceans. On Titan, the surface 
warms quickly relative to a season, and one 
therefore might expect the convective up- 
draft to follow the latitude of maximum 
solar forcing. However, the atmosphere 
could not quickly respond to such warming 
since the radiative cooli ng time exceeds th e 
length of a Titan year (IFlasar et al.lll98ll ). 
For reasonable estimates of surface heat ca- 



pacity (JTokanoll2005l ) , the radiative cooling 
time of the surface is perhaps a month; it 
is clear that the atmosphere provides the 
dominant source of heat capacity in the cli- 
mate system. 

Dry (i.e. neglecting methane condensa- 
tion and evaporation) GCM simulations 
suggest the large heat capacity of the at- 
mosphere does not significantly attenu- 
ate the seasonal cycle of the Hadley cell; 
the convective updraft of the Hadley cell 
travels from one pole to the other each 
Titan year, roughly in phase with the 
seaso nal cycle of solar radiation (ITokano 
20051 ) . However, increasing the amount 
of methane evaporation at the surface in- 
troduces a time lag of the response of 
the Hadley circulatio n to seasonal forcing 
flMitchell et alJl2006[ l. Moist, precipitat- 
ing convection penetrates deeper into the 
cold, dense atmosphere than dry convec- 
tion; one can think of the increase in con- 
vective depth arising from a reduction in 
the vertical temperature gradient resulting 



from moist convection versus dry convec- 
tion (the lapse rate). By doing so, moist 
convection communicates surface warming 
into the large heat capacity reservoir in the 
upper troposphere that is not accessed by 
dry convection alone. The Hadley cell effi- 
ciently communicates convective warming 
throughout the latitudinal domain, which 
increases the thermal inertia of the sys- 
tem. The net result is an increase in the 
mass of the atmosphere through which the 
warming of the subsolar surface is com- 
municated, which introduces a time lag in 
the response of the Hadley cell to seasonal 
forcing. 

Figure [3] shows the predicted vertical 
pressure velocity with the scale truncated 
to highlight regions of large-scale updraft 
(colored contours) as a function of lati- 
tude and date over the last half Titan year 
of three model simulations^!; the insolation 
pattern is overplotted for reference (black 
contours) . The position of the most intense 
updraft features in Figure [3] relative to the 
solar forcing gives a measure of the time 
lag of the Hadley cell. The "dr y" lim it 



resem bles other dry simulations flTokano 



200 5]), and very litt le time lag is predicted 



( ILorenz et al.ll2008l ). In the "moist" case, 
in which deep, precipitating convection is 
most favorable, the main convective up- 
draft is in quadrature with the seasonal 
forcing. The "intermediate" case, in which 
methane evaporation was restricted to pro- 
duce ~50% relative humidities near the 
surface, has sporadic precipitating convec- 
tion, and as a result has a smaller time lag. 
Atmospheric angular momentum (with 



This model is sim ilar to our earlier one 
( Mitchell et al.1 1200(a) except a m ore realis tic ra- 



the minimum subtracted) and the induced 
spin rate drift (for an assumed 70 km ice 
I shell with the average drift subtracted) 
over one Titan year for our three simulated 
cases are shown in Figures [1] and [2, respec- 
tively. The time lag of the Hadley cell up- 
draft evident in Figure [3] translates directly 
into a time lag in these quantities. For 
comparison, the measured spin rate drift 
and the ac celeration are shown as a solid, 



bold line (ILorenz et al.l 120081 ; IStiles et al. 
20081 ). The "intermediate" case fits the 
data very well, indicating methane ther- 
modynamics plays a role in establishing 
the observed time lag in the measured spin 
rate drift. While our dry model repro- 
duces the change in atmospheric angular 
momentum in the preyious model estimate 
(JTokano fc Neubauerll2005l ). the change is 
somewhat smaller in the intermediate and 
moist cases. Our model predicts the spin 
rate drift should currently be changing 
directions; Cassini RADAR observations 
through 2010 should unambiguously deter- 
mine whether our prediction holds. 



4. Thermodynamic efficiency of Hadley 
cell oscillations: Comparison with 
Earth 

Work must be done to change the spin 
rate of the ice shell. The torque exerted 
by Hadley cell oscillations is driven by sea- 
sonal changes in insolation. The thermo- 
dynamic efficiency of this process can be 
estimated given the change in rotational 
kinetic energy and the integrated sunlight 
over a quarter Titan season. 

The kinetic energy of the surface shell is 



T, 



-LW 



diation scheme was implemented (Mitchell 2008) 



(8) 



The change in kinetic energy resulting from 
a change in surface angular momentum, 
AL = I S AQ, is approximately 



AT, 



ALQ . 



(9) 



For the numbers listed previously, AT S ~ 
10 20 Joules. The time to go from maximum 
to minimum spin is roughly a single season, 
or ~7 years; in this time, Titan receives 7 
xlO 22 J of radiant energy from the Sun. 
The thermodynamic efficiency of seasonal 
angular momentu m exchange is e ~ 0.1%. 
Sohl et al.l (120031 ) suggest tidal dissipation 
in Titan could be 10 10 W which over a Ti- 
tan season amounts to ~ 10 18 Joules, two 
orders of magnitude smaller than the sur- 
face kinetic energy change. 

Earth experiences seasonal changes in 
length of day of a millisecond due to at- 
mospheric torques. Earth's moment of in- 
ertia about the spin axis is 8 xlO 37 kg m 2 , 
and the kinetic energy associated with the 
spin ra t e cha nge is ~ 5 x 10 21 Joules (see 
Rosen! (119931 ) for a review of the topic). 
Earth receives ~ 10 24 Joules of radiant 
energy from the Sun in a season, imply- 
ing an efficiency of e ~ 0.5%, the same 
magnitude as that of Titan. Titan re- 
ceives l/15th of the solar energy the Earth 
receives over the time of their respective 
seasons (7 years for Titan, 3 months for 
Earth). Taking the change in AAM as a 
measured quantity, the difference in solar 
energy nearly cancels the factor of 16 due 
to the difference in rotation rate in Equa- 
tion [9] when estimating the efficiency. The 
change in Titan's relative AAM, L max , is 
~l/5th of Earth's, which accounts for Ti- 
tan's lower efficiency. However, since Ti- 
tan's stratosphere either absorbs or reflects 
more than half of the incoming solar radi- 



ation (IMcKay et al.lll99ll ). Titan's tropo- 
spheric Hadley cell is driven by a reduced 
level of solar radiation relative to the top- 
of-atmosphere. If this is taken into ac- 
count, the efficiency increases to near the 
level of Earth. 

Rosenl (119931 ) attributes Earth's sea- 



sonal changes in length-of-day to variabil- 
ity in the subtropical jets, which mark 
the poleward flanks of Earth's angular- 
momentum-transporting Hadley cell. Fig- 
ure H] compares changes in relative AAM 
on Earth fro m NCEP/NCAR re analysis 
data for 2007 flKalnav et al.l Il996h to the 
our "intermediate" Titan simulation; the 
top two panels display the vertical mass- 
weighted average of ua cos(^), which is the 
specific angular momentum of the atmo- 
sphere relative to rotation (measured pos- 
itive in the direction of rotation), and the 
bottom two panels show the correspond- 
ing relative AAM. Earth's subtropical jets 
contribute bands of positive AAM between 
30-40 N/S latitude in Figure 11(a). The 
Northern hemisphere jet has stronger sea- 
sonality, with a minimum during North- 
ern hemisphere summer corresponding to a 
global minimum of relative AAM in Figure 
IUc). Titan's subtropical jets are spread 
further poleward because of the increased 
width of the Hadley cell (Figure I3jb)), but 
the pattern with season shows similarity to 
the Earth's. In particular, minima in rela- 
tive AAM occur when either the Northern 
or Southern subtropical jet is at minimum. 

5. Other potential sources of time 
lag 

There are a few physical processes which 
could additionally contribute a time lag, 
and we presently estimate their effects. 



First, a sluggish atmospheric overturning 
circulation might contribute a time lag; 
we assumed in our analytic model that 
the circulation is vigorous and able to re- 
distribute the angular momentum of the 
convective updraft instantaneously. The 
overturning timescale varies by an order 
of magnitude between the "moist" case 
and the "dry" limit. Since our analytic 
model assumes the atmosphere instanta- 
neously acquires the angular momentum of 
the surface at the updraft, the effect of a 
sluggish circulation is likely positively con- 
tributing to the time lag in AAM. Since the 
atmosphere is most sluggish in the "moist" 
case, the effect goes in the same direc- 
tion as the thermal inertia. The weaker 
Hadley circulations in the "intermediate" 
and "moist" cases reduce the amplitude of 
AAM changes relative to the "dry" limit, 
as seen in Figure 1 in the main text; as a 
result, we infer a somewhat shallower ice 
shell depth of ~70 km. 

Second, frictional (Ekman) drag at the 
interfaces of the subsurface ocean could dy- 
namically couple the surface to the deep 
core. Though the drag is only communi- 
cated to a thin layer of depth S ~ (2v/f), 
where v is the viscosity of water and / ~ 
2Q is the coriolis parameter, mass trans- 
port in this layer converges (diverges) to 
form downwelling (upwelling) through the 
depth. The characteristic timescale of this 
overturning is TEk y Hj^fu) - 1 / 2 , where H 



i s the ocean depth (jVaHisll2006l ). ISohl et al. 



( 120031 ) estimate the ocean depth to be 224 
km for an assumed 15 wt.% NH 3 . If we as- 
sume H ~ 200 km and / ~ 2f2, then TEk ~ 
3000 years. Ekman drag in a laminar 
boundary layer cannot contribute a time 
lag of a few years. However, the Richard- 



son number, Ri, in an assumed laminar 
boundary layer may be order unity}. In 
this case the molecular viscosity, u, should 
be replaced with an effective eddy viscos- 
ity, k. We estimate k ~ L d f /y/lti, where 
Ld ~ u/ f is the radius of deformation 
(I Vallisl 120061 ) and u ~ Af2a is the velocity 



difference in the boundary layer induced 
by libration of the crust. For previously 
listed values, k~2-3x 10~ 2 m 2 /s and 
the Ekman spin-down due to the turbulent 
boundary layer is TEk ~ 20 years. This 
rough estimate suggests that if the bound- 
ary layer is turbulent, dynamic coupling of 
the surface and core by Ekman drag could 
introduce a time lag in the drift rate. 

Third, we have implicitly assumed 
torques on a permanent quadrupole mo- 
ment in the crust cannot overcome fric- 
tional torques provided by the atmosphere, 
which may not be the case. The full equa- 
tion of motion (EOM) of the surface crust 
including these terms reads 



C s & — T atm + Tsat + T, 



coup ) 



(10) 



where C s = I s is the polar moment of in- 
ertia of the crust and = 0.1825 deg/yr 2 
is the measured rotational acceleration of 
the surface. The T's represent torques of 
various mechanisms averaged over an or- 
bital period. T atm is the seasonal torque 
of the atmosphere; our analysis has as- 
sumed this term dominates the others, 
C s ~ T atm . T Sa t is the torque of Saturn's 



*We estimate the interior heat flux from the model 



of lSohl et al.l (|2003l ) and calculate the temperature 
difference over the laminar Ekman layer due to 
thermal conduction. We then use this temperature 
difference along with an estimate of the velocity in 
the boundary layer, u ~ Afla, to estimate Ri ~ 
0.7 



gravitational field acting on permanent de- 
formations contributing a quadurpole mo- 
ment to the crustal layer. T coup represents 
gravitational coupling between a perma- 
nent surface quadrupole and a permanent 
quadrupole in the core. Gravity measure- 
ments from the Cassini spacecraft can con- 
strain the quadrupole moments, but none 
have yet been reported in the literature. 
By assuming they do not contribute sig- 
nificantly to Equation [TUl we have implic- 
itly placed limits on the surface and core 
quadrupoles contributing to T Sa t and T coup , 
which we presently quantify. We will ex- 
press limits in terms of the ratio (B—A)/C, 
where C is the polar moment of inertia and 
B and A are the moments orthogonal to 
the spin axis. The subscripts "s" and "c" 
will respectively refer to moments of the 
surface and core. 

Saturn's gravity would apply a torque 
to a permanent quadrupole in the sur- 
face. Neglecting this effect in the EOM 
implies 6 >> Ts a t/C s . The torque is 
proportional to (B — A) s , the quadrupole 
of the surface, and n ~ Q, the orbital 



mean motion flMurray fc Dermottl Il999l ). 
Together with the Cassini measurement of 
rotational acceleration, the constraint im- 
plies {{B—A)/C) s « 1(T 7 . Assuming the 
surface is a dynamically isolated ice shell of 
70 km thickness, (B - A) s « 2 x 10 27 kg 
m 2 . The maximum (B — A) s could, for in- 
stance, be due to an isostatically uncom- 
pensated disk of ice of depth 1 km and 
diameter smaller than 40 km sitting on 
the equatorial surface; a 100 m thick disk 
would need to be smaller than 125 km. 

A surface quadrupole could also grav- 
itationally couple to a quadrupole in the 
core, represented by (B — A) c . As ex- 



pressed in the EOM, we have assumed 
6 >> T coup /C s . We derive a constraint 
on ((B — A)/C) c by making the follow- 
ing assumptions: 1) the quadrupole of the 
surface layer relative to its polar moment 
takes the maximum value estimated in the 
previous paragraph, ((B — A)/C) s ~ 10~ 7 ; 
2) the core remains fixed in the rotating 
frame; and 3) the core radius is 2250 km 
a nd has the compos ition shown in Figure 6 



of ISohl et al.l (120031 ). Under these assump- 
tions, {{B - A)/C) c « 5 x 10- 5 . This 
minimum value is intermediate between 
the moon (which is roughly the size of Ti- 
tan's rock and iron core), ((B — A)/C) m ~ 
2 x 10~ 4 , and the Earth, {{B - A)/C) & ~ 
2 x 10~ 5 . 

6. Conclusions 



We have developed an analytical model 
for the time-dependent atmospheric an- 
gular momentum in an atmosphere with 
an oscillating Hadley cell which agrees 
with a simulation of Titan 's atmosphere 
(ITokano fc Neubauerl 120051 ); while both 
models predict the observed magnitude of 
spin rate drift, neither can produce the 
observed time lag of the seasonal drift. 
Numerical simulations with methane ther- 
modynamics included show an increase in 
the thermal inertia of the system, thus 
introducing a time lag with the seasonal 
forcing. The time lag is sensitive to the 
supply of methane evaporation at the sur- 
face. If methane is allowed to efficiently 
evaporate from the model surface, deep 
precipitating convection persists at all sea- 
sons and warms the high-heat-capacity at- 
mosphere which increases the thermal in- 
ertia of the system. In this "moist" case, 
the convective updraft of the Hadley circu- 



lation is in quadrature with solar forcing, 
which translates to a ~ 7 year time lag. 
If instead surface evaporation is uniformly 
restricted, the numerical model has the 
appropriate amount of time lag, matching 
the observed lag in the spin rate drift (see 
Figure [2]). Our model predicts the spin 
rate drift has changed directions; current 
Cassini RADAR observations of Titan may 
show this trend. 

In conclusion, the observed time lag of 
Titan's drift in rotation rate is consistent 
with the lag introduced by thermodynamic- 
dynamic coupling resulting from the methane 
cycle in the atmosphere of Titan. We have 
shown the change in total atmospheric an- 
gular momentum is primarily controlled 
by the latitude of the Hadley cell updraft, 
and this updraft marks the location of the 
most persistent clouds features. The re- 
cent disapp earance of pers i stent south po- 
lar clouds (jSchaller et al.l 120061 ) suggests 
the Hadley cell updraft is moving equator- 
ward in response to the changing season. 
The resulting transfer of angular momen- 
tum to the atmosphere from the surface 
should currently be decelerating the spin 
rate, as shown in Figure [21 Other sources 
of torque might measurably contribute to 
the observed spin rate drift; however bet- 
ter knowledge of Titan's gravitational mo- 
ments and interior structure is required to 
estimate their effects. 
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Fig. 1. — Relative atmospheric angular momentum (AAM) over one Titan season for our 
analytic and numerical models. See text for descriptions. 
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Fig. 2. — Induced spin rate drift over one Titan season assuming a uniform, dynamically 
isolated ice I shell of 70 km thickness. Lines are defined as in Figure [U 
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Fig. 3. — Updraft velocity at the lowest model layer in mm/s (colored contours) and solar 
forcing (black contours) during the Cassini epoch for our numerical simulations. See text for 
definitions of each case. 
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Fig. 4. — Vertical (mass-weighted) average of specific relative atmospheric angular momen- 
tum of Earth (a) and Titan (b) as a function of latitude over one full seasonal cycle. N/SSE - 
Northern/ Southern Spring Equinox; N/SSS - Northern/ Southern Summer Solstice. Global 
relative atmospheric angular momentum of Earth (c ) and Titan (d). Ea rth diagnostics per- 
formed on NCEP/NCAR reanalysis data for 2007 ( Kalnay et al.lll996l ). Titan diagnostics 
performed on the "intermediate" simulation. 
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